Optimal thermoelectric figure of merit of a molecular junction 
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We show that a molecular junction can give large values of the thermoelectric figure of merit ZT, 
and so could be used as a solid state energy conversion device that operates close to the Carnot 
efficiency. The mechanism is similar to the Mahan-Sofo model for bulk thermoelectrics — the Lorenz 
number goes to zero violating the Wiedemann- Franz law while the thermopower remains non-zero. 
The molecular state through which charge is transported must be weakly coupled to the leads, and 
the energy level of the state must be of order ksT away from the Fermi energy of the leads. In 
practice, the figure of merit is limited by the phonon thermal conductance; we show that the largest 
possible ZT ~ (G5'h)"'''^ 
conductance quantum. 



where G^^ is the phonon thermal conductance divided by the thermal 



PACS numbers: 73.23.-b, 73.23.Hk, 72.15.Jf, 85.65.+h 
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Electrical transport in mesoscopic and nanoscale sys- 
tems has been a major focus of recent work in physics, 
chemistry, and materials science. This interest has been 
motivated both by the potential for applications in elec- 
tronics and by the opportunity to study novel trans- 
port phenomena such as conductance quantization^ and 
Coulomb blockade. Compared to electrical transport, 
thermal transport has been much less intensively stud- 
ied, despite a number of fascinating results including the 
observation in 2000 of the quantum of thermal conduc- 
tance^' Thermal and thermoelectric transport are also 
relevant to technological questions of great importance, 
such as the construction of solid-state energy-conversion 
devices. 

The efficiency of such devices is determined by the 
thermoelectric figure of merits^ ZT, defined as 



ZT = 



GS^ 



GfjT- 



Here G is the charge conductance, S the thermopower, 
Gj'j the electron thermal conductance, and G^^ the 
phonon thermal conductance. As ZT — > oo, the device 
attains the Carnot efficiency. The figure of merit of the 
semiconductors typically used for thermoelectric appli- 
cations is approximately 1; ZT values of 2 or 3 would 
lead to a considerable increase in commercial utility. In- 
creased interest in ZT has given additional motivation 
to work on the thermal and thermoelectric properties of 
unusual materials: nanowires^, nanostructured mate- 
rials^,^ strongly correlated materials2.i^iiS, and molecular 
junctions"'^^'^^. Of these, molecular junctions are partic- 
ularly promising, as it is hoped that the phonon thermal 
conductance will be small due to a density of states mis- 
match between the phonons in the leads and the phonon 
modes of the molecule. Improvements of standard den- 
sity functional methods are needed to give unambigu- 
ous predictions for the transport properties of molecular 
systemsi^iii^, especially when correlations between elec- 
trons are significant. It is therefore valuable to have a 



thorough understanding of transport in simplified mod- 
els, such as the Anderson model (reviewed below), that 
describe intramolecular correlations, both as a check on 
the accuracy of new ab initio methods and to guide both 
experiment and numerics to where high ZT may occur. 

In this paper we find the figure of merit of a molecu- 
lar junction. We start from the Anderson model for the 
molecule and leads, which is valid under the assumption 
that a single molecular level dominates transport. If the 
lifetime r of an electron in the molecular state satisfies 
T < h/{kBT), it is generally truei^>i^ that ZT < 1. We 
concentrate therefore on the regime t ^ h/ [kBT). Phys- 
ically, this corresponds to a wavefunction for the level 
that is strongly localized at the center of the molecule, 
so that the state is weakly coupled to the leads. In ref- 
erence (fl7h it was shown that for the Anderson model in 
this limit of weak coupling to the leads, with strong in- 
teractions, and ignoring the phonon contribution to the 
thermal conductance, the figure of merit becomes infi- 
nite. Here we find that the phenomenon is surprisingly 
general — it is also true in the limit of weak electron- 
electron interactions. We clarify the mechanism in both 
cases by comparing with the work of Mahan and Sofc^* 
on bulk thermoelectrics. (The mechanism, a violation 
of the Wiedemann- Franz law, has previously been stud- 
ied in quantum dots in different parameter regimes in 
references (HSiOill).) Also, by carrying the calculation 
for the electron-only ZT to higher order in the coupling 
to the leads, we find the optimal value of ZT once the 
phonon contribution is included. Finally, we suggest a 
candidate molecule where high ZT might be observed. 

The Anderson model is defined by the Hamiltonian 

H = -t ^ (c^j+i_o.Cj_CT + he) 

J>0j<-1,CT 

~ {4(c^,<tCo,^ + he) + t'L{c^o,aC-i,a + he)} (1) 

(T 

The label a represents the spin of the electrons; j is the 
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site index, with the molecule represented by the site at 
j = 0. The amplitude to hop from the molecule to the 
right (left) lead is given by t'^^j^y and the leads are rep- 
resented by the sites at j > and j < 0. The energy to 
place a single electron on the molecule is Ed, which can 
be interpreted as a gate voltage. The electron-electron 
repulsion within the molecule results in an additional en- 
ergy cost of U to add two electrons to the molecular site. 
The hopping within the leads is t, which gives a band 
structure of E{k) — — 2tcosfca, where a is the lattice 
spacing. 

As a first step, we find the transport properties with no 
electronic repulsion ([/ = 0), when the transmission func- 
tion of the impurity can be found by solving the single- 
particle Schrodinger equation. The result is 



i?2/(4i2)) + - r/(2t)) - £d)2 



where E is the energy of the incident electron, ^l{r) = 

2^'i(_R)Aj and T = {Tl + r_R)/2. For electrons close to 
the center of the band, and Tl^h <^ t, we get the familiar 
Lorentzian transmission 
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The hybridization energy, F, determines the mean life- 
time, r, of an electron on the dot, through r = h/T. 
We will assume this form of Ta{E) in the study of the 
non-interacting case that follows. 

Within the Landauer formalism the transmission func- 
tion determines the charge current / and the heat current 
due to the electrons Iq through the relations 



ret 

J^n — 



dE{T^{E)+T^{E)){fl{E)-f%{E)l 

-2t 

dE{E t,){T^{E) + T,{E)){fl{E) - fl{E)). 
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Here — e is the electron charge, and the function /° is 
defined by fl{E) = {^e'^E-t,r.)/{kBTr.) + .^^here 
and Tl are the chemical potential and temperature in 
the left lead (and similarly for /]^). At linear response, 
and for a chemical potential away from the band edges, 
we have 

/ - -^2!^ (^0(7, 5)eV + ^1(7, 5)kB^T) , (2) 
/q' = (J-i(7, 5)eV -f ^2(7, 5)kBm , (3) 



where 7i(^) = T l^^)/ {ksT), 7 = r/{kBT), 5 ^ {sd - 
fi)/(kBT), eV = hl - and AT ^ Tl - Tr. The 
functions can be found from equations ^ and ([3]) by 
setting X — {E — iJ,)/{kBT)\ the result is 



Ml, 6) 



dx 
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{2cos\Y{xl2)f^^ + {x^5Y 




FIG. 1: The Lorenz number (main figure) and the ther- 
mopower (inset) plotted as a function of T/(kBT), with 
{ed — l-i)/ (kBT) = 2 in both cases. Only the electron contribu- 
tion to the thermal conductance is included in the calculation 
of the Lorenz number. 



J-Q, T\ and J-2 can then be expressed as 



I J /tt + w 
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where w = 7 + j(5. The function -0 is the trigamma bmc- 
tion il^iz) = Y.n=oi^ + "-)~^- 

We define the limit of a long-lived molecular state as 
F ^ fcsT, i.e. 7 <C 1. Ignoring the phonon contribution 
to the heat current (to be restored below), the figure of 

merit is given by ZT — {TqT2IT\ — l) . For 7 fixed, 
ZT is largest when 5 ~ ±2.4. If one holds 8 fixed, and 
allows 7 — > 0, the figure of merit diverges as 



ZT = 



1 



0(7"). 



7 (2cosh((5/2))2 

(This follows from expanding JF„ to first order in 7.) As 
shown in Fig. [Tl the thermopower is finite as 7 ^ 0, 
and the divergence in ZT is due to a violation of the 
Wiedemann-Franz law, with the thermal conductance 
vanishing faster than the electrical conductance. 

In practice, the figure of merit of the junction will be 
limited by the phonon thermal conductance. This can be 
parametrized in terms of the fraction G^^ of a thermal 
conductance quantum transmitted by the molecule, i.e. 
Using 



nv^ — nv^ 



th 3 h 



GS' 



ILIR 2fc| ttJ^ 



7 h (2cosh(5/2))2^ 
we include the phonon contribution to get 
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FIG. 2: 1/ZT, without the phonon contribution to the ther- 
mal conductance, is plotted for T/(kBT) and {sd — 
li)/{kBT) = 2. 

Assuming 7^ = 7_r, — 7 and minimising l/ZT with re- 
spect to 7 and S gives, for small G^/j, an optimal value 
of 

ZT~0.51/(G^;^)i/2. 

The figure of merit becomes infinite as G^^ goes to zero, 
as we found before. 

There is a simple physical picture that gives us the cor- 
rect expressions for JFqi ^1, and to order 7°, and thus 
gives 1/ZT = 0-1-0(7''). This approximation is known as 
"sequential tunneling" , and has been used to study trans- 
port in quantum dots and molecule o^^'^'^i^'*'^^ . Suppose 
the molecular level has energy Sd- Then, in this approxi- 
mation, an electron in the leads can tunnel onto the level 
only if it too has energy Sd- Energy is therefore conserved 
throughout the transport process. This allows for a cur- 
rent proportional to TlTr/T. In quantum-mechanical 
tunneling by contrast, the intermediate state is virtual; 
away from resonance, the current is proportional to the 
probability to tunnel onto the level multiplied by the 
probability to tunnel off, and thus to r^ri^. 

In the sequential tunneling limit, each electron passing 
through the junction carries the same amount of heat, 
Q = Ed — fJ-- To understand why the Wiedemann-Franz 
law is violated, it is essential to bear in mind that the 
thermal conductance is defined as the heat current di- 
vided by the temperature difference under the condition 
of zero charge current. The condition / = implies that 
the flux of left-moving electrons is the same as the flux 
of right-moving electrons; since each left or right moving 
electrons carries exactly the same amount of heat Q, the 
thermal current Iq is also zero. This is the same mecha- 
nism that Mahan and Sofo used to model the bulk mate- 
rial with the best thermoelectric figure of meri t^^i^^ . For 
long-lived molecular states, unlike bulk materials (except 
in the atomic limiti^), the theoretical analysis can be 
extended to include electron-electron repulsion, as now 
done. 



We now consider whether, once electron-electron re- 
pulsion in the molecule is added, the electronic figure of 
merit ZT still scales like 1 /7 in the limit of a long-lived 
molecular state. Define the probability of the molecu- 
lar level being empty as ^'(0), of there being an up spin 
as -P(t), a down spin as P{i), and of it being doubly 
occupied as P{2). We consider the currents where the 
molecule meets the left lead. An example of a process 
that contributes to the currents at this point is the tun- 
neling of an up electron with energy Sd from the left lead 
onto the empty level. To find the total currents, we con- 
sider all possible tunneling processes, both onto and off 
of the molecule. For S > 0, ie. Ed > fi, the charge and 
heat currents at the left junction are then given by 

/ = -er-i (P(0)2/L(ed) - P(l)(l - hied)) 

+Pil)fL{ed + U)- 2P(2)(1 - fLied + U))) , 
= T^'keT idPrnfUsd) ~ SPim - hisd)) 
+ {5 + u)P{l)fL{ed + U) 
-{5 + u)2P{2){l-fL{ed + U))), 

where P(l) = P(t) -hP(i), u = U/iksT), and, as before, 
5= {sd- fi)/ikBT). The probabilities F(0), P(t), P{i), 
and P{2) are determined from the steady state condi- 
tion. The final result for the linear response transport 
coefficients is 

^_e^ lLlR e^+"(l-He^(2 + e^+")) 

h 7 (1 + 6-5) (1+6-^+") (1+6-5+" (2 4-6-5))' 

fca (<5(l + 6^(2 + e^"+")) + M(l + 6'^)) 
e (1 + 6-5(2 + 6'5+«)) 

G|[ ^ fclTLTfl u'e^'+" 

T ft 7 (l + 2e*+« + e2-5+'')(l + 2e-5 + e2-5+'')' 

(It is straight-forward to confirm that this result, with 
u = 0, agrees with the exact result for the non- interacting 
case in the limit of jl,fi. ^ 1, if one makes the identifica- 
tion Ti/^ji-) = ^L{R)/f^- This confirms that only sequential 
tunneling is important in the limit of a long-lived state.) 

The figure of merit in the sequential tunneling approx- 
imation (again ignoring phonons) is given by 

_ 6-^ (.5(1 + 2e^" + 6^^+") + u{l + 6^))^ 
u2 (i + e*)(l + 6-5+") 

(A n approximate version of this formula was derived in 
[13.) l/ZT is plotted in Fig. Has a function of u for (5 = 2 
and 7 — > 0. At u = 0, the non-interacting case, ZT is di- 
vergent, as before. For u 3> 1, 1/ZT goes to zero like e~". 
This can be understood within the picture presented ear- 
lier. Even though, in contrast to the non-interacting case 
there are now two levels (of energy Sd and Sd + U) on the 
molecule through which transport can occur, the relative 
probability of electrons in the leads having the corre- 
sponding energies goes as 6^". As U becomes large there 
is thus only an exponentially small probability that there 
will be an electron in the leads with energy Sd + U. In 
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this case, each electron that passes through the junction 
carries heat Sj — fJ- and the situation is similar to the non- 
interacting case where transport takes place through just 
one level on the molecule. Thus, the thermal conductance 
in this case is also zero. (For convenience, let us consider 
the symmetric case, with 7l = 7fl = 7-) Since the se- 
quential tunneling model is correct only to first order in 
7, we can conclude that the thermal conductance van- 
ishes to first order in 7, and so must be of order 7^. For 
i5 ~ 2.4 and m » 1, we have GS^ ~ ^OAkg/h, and so 
GS^ scales like 7 as in the non-interacting case. Adding 
in once more the phonon contribution to the heat current, 
we can again conclude that ZT will scale like (G^^)"^/^. 
In the limit of u ^ 00, the Green function method^liS^ 
can be used to find the transport coefficients to higher 
order in 7, and so the coefficient in the scaling of ZT. 
The result is that, including the phonon contribution as 
before, the optimal value of ZT is 0.42(G'^^)-i/2. 

We have shown that transport through a single molec- 
ular level can lead to large values of the thermoelec- 
tric figure of merit ZT. Within the Anderson model 
for the level, this occurs if the molecular state is long- 
lived (i.e., for the case of a symmetric junction, 7 — > 0) 
and the energy scale U of the electron-electron interac- 
tions within the molecule satisfies either U <^ ksT or 
U ^ ksT. The increase in ZT is due to a violation 
of the Wiedemann- Franz law — the charge conductance 
scales like 7, whereas the electron thermal conductance 
divided by the temperature scales like 7^. While sim- 
ilarly large figures of merit may be achieved for opti- 



mal band structures in bulk materials*^, the advantage 
of molecular systems is that they may have only a small 
phonon contribution to the thermal conductance, so that 
their maximum achieved value of ZT once phonons are 
included may be higher than in bulk materials. 

There are some practical challenges in achieving dra- 
matically enhanced ZT in a real device by the mecha- 
nism presented here. The wavefunction of the level must 
be localized at the center of the molecule to give a small 
overlap with the states in the leads, and so a long lifetime; 
in addition, the energy of the level must lie a distance of 
about 2kBT from the Fermi level of the leads. These 
requirements may be satisfied by the Co(tpy-(CH2)5- 
SH)2 molecules studied by Park et a P^i'^° . The electronic 
transport in these systems is through a cobalt atom that 
sits at the center of the molecule; the alkyl chains sepa- 
rate the cobalt from the gold leads, creating a state lo- 
calized away from the leads. It would be of great interest 
to see if indeed the Wiedemann-Franz law is violated in 
these molecules, as implied by the results in this paper. 
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